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By making use of Numerical Stochastic Perturbation Theory (NSPT) we can compute renorma- 
lization constants for Lattice QCD to high orders, e.g. three or four loops for quark bilinears. 
We report on the status of our computations, which provide several results for Wilson quarks and 
in particular (values and/or ratios of) Zy, Za, Zs, Zp. Results are given for various number of 
flavors {nf = 0,2,3,4). While we recall the care which is due for the computation of quantities 
for which an anomalous dimension is in place, we point out that our computational framework 
is well suited to a variety of other calculations and we briefly discuss the application of NSPT to 
other regularizations (in particular the Clover action). 
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1. Introduction 

Numerical Stochastic Perturbation Theory (NSPT) is a numerical tool which enables us to 
perform high order computations in Lattice Perturbation Theory (for a comprehensive introduction 
see [[T]]). In recent years our group undertook the high order computation of Lattice QCD renor- 
malization constants. We have till now mainly focused on Wilson quarks bilinears [Q, ||], a task 
which is by now almost completed ([Q] will contain a comprehensive report). Here we report on 
the three (or even four) loop calculation of finite renormalization constants. The overall strategy of 
the computation will be sketched in order to make this communication sufficiently self consistent. 



2. Computational setup 

Our computations are for Wilson gauge action and (unimproved) Wilson fermions. We can rely 
on a collection of NSPT configurations for various number of flavors: = (three loops, both 
on 32'* and on 16^ lattices; the latter is mainly used to check finite size effects); = 2 (on a 32^ 
lattice, to three and even four loops); «/ = 3 and = 4 have at the moment only a reduced sample 
of configurations, which enables us to quote preliminary results (i.e. only for certain quantities). 

The scheme we adhere to is what is known as RI'-MOM, which became very popular after its 
non-perturbative implementation [||]. It is important to point out that a comprehensive computation 
of the quark bilinears anomalous dimensions for this scheme is available in the literature ||^. 

In order to implement our program we start from the computation of quark bilinears between 
external quark states at fixed (off-shell) momentum p 

J dx{p\ \j/{x)r\i/ix) \p) = Gr{pa), (2.1) 

where F stands for any of the 16 Dirac matrices singling out the S,V,P,A,T currents. We notice that 
on the lattice the dependence is on pa. These quantities are gauge-dependent and our computations 
are in the Landau gauge, which is easy to fix on the lattice and does not require to discuss the 
gauge parameter renormalization. Next step is trading the Gr{pa) for the amputated function, in 
our notation Fr (pa) 

Gr{pa) Tripa) = S^^ (pa) Gr{pa) S^^ (pa). (2.2) 
The Fripa) are finally projected on the tree-level structure by a suitable operator P^r 

Or{pa) = TT{Po^rr{pa)). (2.3) 
Renormalization conditions are given in terms of the Or{pa) according to the master formula 



Zoril^a,g{a)) Z \na,g{a)) Or{pa) 



1. (2.4) 



The dependence of the Z's on the scale jj. is via the dimensionless quantity jxa, while the depen- 
dence on g{a) explicitly displays that our results will be expansions in the lattice coupling. The 
quark field renormalization constant (Z^) entering the above formula is defined by 

. 1 Tr(^ S-\pa)) 



Zq{}la,g{a)) 



12 p 



,2 



(2.5) 
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In order to get a mass-independent renormalization scheme, one gives the renormahzation con- 
ditions in the massless Umit. In Perturbation Theory this imphes the knowledge of the relevant 
counterterms, i.e. the Wilson fermions critical mass. One and two loops results are necessary in 
order to implement our three loops computations and are known from the literature As a side 
remark we notice that third (and fourth) loops have been computed by us as a (necessary) byproduct 
of the current computations. We stress that in our (perturbative) framework reaching the massless 
limit does not require any extrapolation procedure. 

At any loop order L, the expected form of the U'^ coefficient of a renormalization constant is 



i.e. we have to look for a finite number, a divergent part which is a function of anomalous dimen- 
sions 7's and irrelevant pieces, which we expect to be compliant to hypercubic symmetry. The 
anomalous dimensions which are needed are to be taken from the literature and subtracted. After 
such a subtraction a hypercubic symmetric Taylor expansion enables us to fit the irrelevant pieces 
given by F{p). It is instructive to write down an example: for the scalar current our master formula 
Eq. reads (in Landau gauge the quark field has zero anomalous dimension at one loop) 



Ms the quantity actually measured from our simulations; after subtracting the log one should 
now proceed to fit the irrelevant pieces in order to get the constant zi^^. This subtraction procedure, 
as we now briefly comment, is actually affected by finite size effects which are to be assessed. 

It is very important to stress that RI'-MOM is an infinite volume scheme, while we are neces- 
sarily computing on finite volumes. This results in a taming of the log, due to pL effects (notice 
that pL = pN, where N is the number of lattice sites). This has been actually verified by comparing 
16^ and 32^ results. The correct one loop result is recovered once one computes the finite size cor- 
rections to the log on a finite volume in the continuum: for a preliminary discussion the interested 
reader is referred to 

3. Wilson quark bilinears to three (and even four) loops 

Having pointed out the delicate issue which is in place in the case of anomalous dimensions, 
we now proceed to discuss the finite cases, for which we do akeady have definite results. 

3.1 Ratios are always safe 

Results to three loops for the finite ratios of renormalization constants Zp/Zs and Zv/Za are 
reported in Table 1 . The quark field renormalization constant present in Eq. cancels out in the 
ratios, together with the divergences that affect Zp and Zj separately. In Fig. 1 we present (one loop) 
examples of the extrapolation procedure which is entailed in fitting hypercubic symmetric Taylor 
expansions: powers of p = pa are the irrelevant contributions to be discarded in the continuum 
limit, in which the constant (intercept atp = 0) is singled out. Notice that for Zy /Za one can inspect 
families of curves: for given values of p^ they correspond to different lengths of momentum p in 



L 



ZL = CL + Y.dimogipy+Fip) {p = pa) , 



(2.6) 




(2.7) 
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Figure 1: Computation to one loop of finite ratios of renormalization constants: Zp/Zs (left) and Z^jZy 
(right). For the latter different curves correspond to different lengths of the component of momentum p 
along the relevant direction (see text). Solid lines join the data points which have been taken into account in 
the fit. 



the relevant direction \i {e.g. for Y7^¥)^ they correspond to different irrelevant contributions 
(for a short discussion see [^). 

Many checks have been performed on the expansions. Finite size effects are well under con- 
trol, as checked by comparing results on 16^ and 32^ lattices. Notice also that one can compute 
(for example) both Z^/Zy and Zy jZ^. Different ratios come from different (although correlated) 
combinations of data (this is due to the order by order nature of the computation): to a very good 
precision the series obtained are inverse of each other. 



ZpjZs 








Hf 


0(/3-i) 









- 0.487(1) 


- 1.50(1) 


- 5.72(3) 


2 


- 0.487(1) 


- 1.46(1) 


- 5.35(3) 


3 


- 0.487(1) 


- 1.43(1) 


- 5.13(3) 


4 


-0.487(1) 


- 1.40(1) 


- 4.86(3) 


Zv/Za 








"/ 











- 0.244(1) 


- 0.780(5) 


- 3.02(2) 


2 


- 0.244(1) 


- 0.759(5) 


- 2.83(2) 


3 


-0.244(1) 


- 0.744(6) 


- 2.72(2) 


4 


-0.244(1) 


- 0.732(6) 


- 2.57(2) 



Table 1: The finite ratios of renormalization constants Zp/Zs and Zy /Za for various number of flavour n/. 
The four loops results which are available for n/ = 2 are not quoted. 
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One can also check that the rif dependence is (within errors) compatible with the expected one 
(linear at two loop, quadratic at three loop): see Fig. 2. 

To proceed to four loops one needs to plug in the three loops critical mass counterterm. Since 
we know it to a sufficiently good accuracy for the case «y = 2, in this case we also have preliminary 
four loop results, which we do not quote in this communication. 

3.2 Zy and Za 

One loop examples of computations of Zy and Za are depicted in Fig. 3. Zy and Za are finite by 
themselves. In our master formula Eq. they are interlaced with log's coming from the quark 
field renormalization constant. The latter can be eliminated either directly from measurements 
of the propagator or by taking ratios with the conserved vector current. Both procedures return 
consistent results, which are summarized in Table 2. 




0.2 0.4 0.6 0.8 1 '0 0.2 0.4 0.6 0.8 1 

ipaf Ipnf 



Figure 3: Computation to one loop of finite renormalization constants: Zy (left) and (right). As in 
Fig. 1, different curves correspond to different lengths of the component of momentum p along the relevant 
direction (see text) and solid lines join the data points which have been taken into account in the fit. 
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Zv 








«/ 


0(j8-i) 


0(/3-2) 


0(/3-3) 





- 0.487(1) 


- 1.50(1) 


- 5.72(3) 


2 


-0.487(1) 


- 1.46(1) 


- 5.35(3) 


Za 








Hf 






0(/3-3) 





- 0.244(1) 


- 0.780(5) 


- 3.02(2) 


2 


- 0.244(1) 


- 0.759(5) 


- 2.83(2) 



Table 2: The finite renormalization constants Zy and for various number of flavor nf. The four loops 
results which are available for «/ = 2 are not quoted. 



Once again, also 4 loops results are available for rif = 2. Fig. 4 displays resummations of 
Za and Zy at jS = 5.8. In this resummations we do take into account also four loops results. In 
the spirit of Boosted Perturbation Theory (BPT) [||], we show different resummations for different 
definitions of the coupling (P is the basic plaquette): 



1 



log(P) X2 = 



(3.1) 



Our («/ = 2, j8 = 5.8) results are Za = 0.79(1) and Zy = 0.70(1) (one can compare to [g]). Notice 
how it is impossible to assess the consistency of a 1 loop BPT result without knowing higher orders 
in the expansion. 
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Figure 4: Resummations of (left) and Zy (right) for = 2 at j3 = 5.8 to one (circles), two (squares), 
three (diamonds) and four (triangles) loops. We show resummations for different couplings: on the x-axis, 
the (different) values of the different couplings. From the left: xo,xi, X2 (xq is j3 " ' , see text for the definitions 
of the other coupUngs). 



6 



Renormalization constants for Lattice QCD: new results from Numerical Stochastic Perturbation Theory 
Francesco Di Renzo 



4. Work in progress 

An obvious step forward is the computation of quark bilinears in the case of other actions. A 
first possibility is to change the gauge action: we will take into account the Symanzick improved 
version. Some comments can be made for the fermionic regularization dictated by the Clover 
action, for which the knowledge of csw is required. The one loop result is available, which suffices 
for two loops computations. Still, an accurate effective parametrization of csw both for «j = and 
= 2 is available, which entails the one loop result [[l0|]. For example, in the case = 2 

1 - 0.454gg - 0. 175gg + 0.012g^ + 0.045gf, 
'''' = l-0.720g2 ^^-'^ 

Our strategy is of course to compute the second loop coefficient for csw- Estimates of three loops 



coefficients can in the meantime be obtained by extracting from Eq. ( pt. ID (and its quenched coun- 
terpart) an estimate for the two loops csw (the rational form of Eq. (11) is well suited to this 
purpose). 
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